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ABSTRACT 

We formulate a vibrationally resolved kinetics for molecular hydrogen and its cation in the primor- 
dial Universe chemistry. Formation, destruction and relaxation processes for each vibrational level 
are studied and included as chemical pathways of the present model. The fractional abundance of 
each vibrational level as a function of the rcdshift is given: a strong deviation from the Boltzmann 
distribution is found at low z. A discussion of the results is provided, also evaluating the effects of 
relaxation processes on the level populations. Analytical fits for some LTE rate coefficients are given 
in the Appendix. 

Subject headings: cosmology: early Universe — methods: numerical — physical data and processes 



1. INTRODUCTION 

The primordial gas after recombination was composed 
of neutral hydrogen and helium atoms with traces of 
electrons, deuterium and lithium, exposed to the radi- 
ation field of the cosmic background radiation (hereafter 
CBR), and in adiabatic expansion. In these conditions, 
only a limited number of gas-phase reactions was possi- 
ble. The first molecules started to form only ten thou- 
sand years after the Big Bang, at a resdhift z « 2000, 
reaching a freeze-out abundance at z « 100: about 10^^ 
for H2, followed by HD, HJ. IIcII+ and several less abun- 
dant species (see Lepp et al. 2002 for a review). 

Despite the simple chemical composition, the kinetics 
of elementary processes in the early Universe is very com- 
plex, due to the presence of a large number of quantum 
states for atoms and molecules involved in the network 
processes, including the interaction with the CBR. De- 
tailed calculations of molecule formation in the early Uni- 
verse have been recently performed by Hirata & Padman- 
abhan (2006, hereafter HP06), Puy & Signore (2007), 
Glover & Jappsen (2007), Schleicher et al. (2008, here- 
after SOS), and Vonlanthen et al. (2009, hereafter V09). 
In addition to these studies, comprehensive collections of 
analytic fits of chemical reaction rates were also given by 
Anninos & Norman (1996), Stancil et al. (1998), Galh 
& Palla (1998, hereafter GP98) and, more recently, by 
Glover & Abel (2008). 

In studies of primordial chemistry, the rovibrational 
levels manifold of the molecules has been usually ignored, 
mainly because of the lack of state-to-state resolved cross 
sections (or, equivalently, reaction rates) of the relevant 
chemical processes. Early estimates of the vibrational 
distribution function of H2 attempted by means of physi- 
cal order-of-magnitude arguments (Khersonskii 1982) led 
to the conclusion that the vibrationally excited states of 
H2 are rapidly deactivated through radiative decay. This 
is because the H-H2 collision time is much longer than 
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the radiative lifetime of vibrational levels, so that the 
population of the vibrational levels is essentially deter- 
mined by the Einstein coefficients and the probability of 
formation of vibrationally excited H2 molecules. HP06 
were the first to compute in detail the vibrational level 
populations of HJ. However, they did not consider col- 
lisional processes inducing vibrational relaxation. 

Because of these simplifications, chemical reaction 
rates arc usually evaluated under the hypothesis of lo- 
cal thermal equilibrium (hereafter LTE). However, non- 
equilibrium distribution function can significantly affect 
the results of such calculations, as shown explicitly by 
Capitelli et al. (2007) for the dissociative attachment 
process of H2. This evidence suggests that the eval- 
uation of the exact vibrational distribution function is 
necessary to model the plasma kinetics, as pointed out 
by Lepp et. (2002). The relatively simple internal struc- 
ture of atomic and molecular hydrogen, and its molecular 
cation HJ, makes possible such a state-to-state approach 
at least for these species. 

In the present work, we follow for the first time the 
chemical composition of the primordial gas, determin- 
ing the population of vibrational levels of H2 and HJ, 
adopting a detailed network of formation and destruc- 
tion channels and relaxation processes, such as radiative 
and atom/ion- molecule collisions. We ignore any con- 
tribution due to subsequent astrophysical processes that 
affect the chemical composition of the gas such as primor- 
dial massive supernovae, which, besides providing heavy 
elements, may also contribute to the formation of H2 in 
the early Universe (Cherchneff & Lilly 2008). 

The paper is organized as follows: in Section 2 we dis- 
cuss the relevant vibrationally resolved reactions and we 
compute state-to-state reaction rate coefficients; in Sec- 
tion 3 we describe our chemical network; in Section 4 we 
present our results for the evolution of the vibrational 
level populations of H2 and H J ; finally, in Section 5 we 
summarize our conclusions. In the Appendix, we give a 
list of fitting formulae for specific rate coefficients under 
LTE conditions. 

2. VIBRATIONALLY RESOLVED REACTION RATES FOR 
H2 AND 

In a hydrogen plasma, H2 and HJ form by various col- 
lision processes, generally in vibrationally excited states. 
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TABLE 1 

ViBRATIONALLY RESOLVED REACTIONS INCLUDED IN THIS WORK 



Chemical process 
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H2(i;) + hi/ ^ H + H 
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H2(i;) + e" ->• H2(v') ^e' + hv 
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hJ(i;) ^- H^(i;') + hv 
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HcH+ +H H+(i') + Hc 
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References: (1) Esposito et al. (1999, 2001); (2) Ci'zek et al. (1998); 
(3) Krstic et al. (1999, 2002), Krstic (2002); (4) Krstic (2003); 
(5) Flannery et al. (1977); (6) Dunn (1968); (7) Allison & Dal- 
garno (1969, 1970), Dalgarno & Stephens (1970); (8) Takagi (2002); 
(9) Celiberto et al. (2001); (10) Wolniewicz et al. (1998); (11) Posen 
et al. (1983); (12) Linder et al. (1995) (see S08). 

For example, in the case of HJ(w), the electron impact 
ionization process produces a known population of vibra- 
tional states given by the corresponding Franck-Condon 
factors for the transitions. For other formation pro- 
cesses, Hj(t!) is created with different populations of vi- 
brational levels. Additionally, inelastic collisions with 
other plasma constituents strongly affect the evolution 
and thermalization of vibrational levels of H2 and H^. 
This fact underlines the importance of such a vibra- 
tionally resolved kinetics, since a priori assumptions and 
approximations could fail in evaluating the vibrational 
level distribution. 

Typically, coUisional processes in a low-temperature 
plasma involve all molecular degrees of freedom (elec- 
tronic, vibrational and rotational); however, in order to 
model the kinetics of vibrational levels of H2 and HJ. 
rotational equilibrium has been assumed. 

For processes not involving H2 and H^, wc have 
adopted the chemical network recently developed by SOS, 
updating the value of the rate coefhcient for the H~ pho- 
todetachment to include the contribution of non-thermal 
photons due to cosmological H and He recombination 
(see HP06). The corresponding rate coefficient for this 
process is given in Table |31 For H2 and H J , we com- 
puted state-to-state rate coefficients from available cross 
sections, including all vibrational levels of their funda- 
mental electronic states (15 and 19 states, respectively). 

A list of the chemical processes included in our calcula- 
tion is shown in Table [T] Each process will be described 
separately in the following subsections. 

2.1. Associative detachment 

The importance of the reaction of associative detach- 
ment for H~ , 

H + H- ^H2(w)+e-, 
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Fig. 1. — H-H~ associative detachment rate coefficient. Trian- 
gles, total rate coefficient obtained from the cross sections of Cfzek 
et al. (1998), according to Eq. JTJ; dashed curve, fit by GP98 of 
the total rate coefficient from Launay et al. (1991); solid curve, fit 
reported in Table [3] 

was first pointed out by Dalgarno (quoted by 
Pagel 1959). Various studies have been published on the 
vibrationally resolved cross sections and/or rate coeffi- 
cients for this process (Bieniek & Dalgarno 1979, Launay 
et al. 1991), underlying the importance of this channel 
of H2 formation (e.g. Dalgarno 2005; Flower 2000). 

Our calculations are based on the detailed data by 
Cfzek et al. (1998). These authors developed an im- 
proved nonlocal resonance model for the description of 
the nuclear dynamics of the collision complex, and 
calculated both associative detachment and dissociative 
attachment cross sections. They provide data for dif- 
ferent values of the angular momentum J of the H^+H 
collision^. For each initial angular momentum, the H2 
molecule can be produced in states with the rotational 
quantum number equal to J + 1 or J — 1 due to selection 
rules; states with the same J are not allowed by sym- 
metry conditions on the total wavefunction and higher 
changes in J are neglected in their model. The total 
rate coefficient, obtained as a sum of the vibrationally 
resolved rate coefficients, 

^tot(Tgas) ~ ^ ^ fc^) (^gas)i (1) 

V 

is shown in Fig. [1] compared with the rate coefficient 
computed by Launay et al. (1991), as fitted by GP98. 
The new fit is given in Table [H Using this rate coeffi- 
cient in a non- vibrationally resolved chemical network of 
the primordial Universe has a negligible impact on the 
H2 abundance. Indeed, this reaction is also one of the 
most important H~ loss channel. Therefore increasing 
the associative detachment rate coefficient cannot lead 
to a production rate of H2 faster than the H~ produc- 
tion rate. For the same reasons, H~ fraction is reduced 
accordingly, by a factor of 4. 

2.2. Photodissociation o/HJ(i;) 
The photodissociation of HJ(i'), 

H+(i)) + /ii^^H + H+, 

^ Data in electronic format are available at 
Ihttp : //utf .mf f . cimi . cz/~ cizek/AD_H2/| 
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Fig. 2. — photodissociation rate coefficient. Solid curve, 
total rate in LTE, summed over all possible initial, from the cross 
sections by Dunn 1968; dashed curve, fit by GP98 based on data 
by Argyros (1974) and Stancil (1994). 

is a typical example of direct process of photodestruction. 
Lebedev et al. (2000, 2003) provided a significant analy- 
sis of bound-free (including also photodissociation) and 
free- free transitions for Hj(i;), limited however to sums 
on the rotational and vibrational quantum numbers. In 
this work, we have adopted the vibrationally resolved 
cross sections calculated by Dunn (1968), evaluated fix- 
ing the rotational quantum number J ~ I. The compar- 
ison between the LTE rate coefhcient obtained using the 
vibrationally resolved data by Dunn and the fit by GP98 
of data by Argyros (1974) for 2500 K< T^ad < 26000 K 
and by Stancil (1994) for 3150 K< T^ad < 25200 K is 
shown in Fig. [2] 

2.3. Radiative association of H and 

The reaction of radiative association of H and H+ , 

H + H+ (v) + hi^, 

has been the subject of several studies (e.g., Ramaker & 
Peek 1976, Shapiro & Kang 1987, Stancil et al. 1993). 
However, no vibrationally resolved data are available in 
literature. For this reason, we have applied the principle 
of detailed balance to the inverse Hj(i;) photodissocia- 
tion reaction, 



„rad.ass. 



7 (2) 

where T^ad is the temperature of the CBR, hvy is the 
thermal threshold for the v-th level. The partition func- 
tions Z of reagents and product are 



Z 



{2TrkTg, 



\3/2 



^.9h(3h for H, 



/l3 



3/2 

%+ 
3/2 



5h+ 



for H^ 



jj+.9H+'9H+^rot(t',rgas) for H+(u), 

where gn, <?h+ and gj^+ are the multiplicities of the nu- 
clear spin contributions for H, and H^, equal to 
5h = 5h+ 2 X (1/2) + 1 = 2 and g^^ = 4; 

and Qtt+ are the electron spin degeneracies, equal to 2 

both for H and under the hypothesis that only their 
ground electronic states contribute to the partition func- 
tion and Ziot(w, Tgas) is the rotational partition function. 




T|K] 

Fig. 3. — H-H+ radiative association rate coefficient. Solid 
curve, total rate coefficient obtained using the detailed balance 
principle from the photodissociation data by Dunn (1968), accord- 
ing to Eq. l^J; dashed curve, fit by GP98 of the total rate coefficient 
from Ramaker & Peek (1976). 

that takes into account the rotational structure of each 
vibrational level of the cation hydrogen molecule, given 
by 

ZUv,T,,s)=J2gje-^''"~''-'^/''^^-, 
J 

where Eyj is the binding energy of the (v, J)-th level 
and Eq the binding energy of the lowest rovibrational 
level. The symbol gj represents the multiplicity of each 
electronic-rotational level. 



^ i 4 



3(2J- 



1) 
1) 



for J even, 
for J odd. 



A comparison between the present LTE calculation and 
the fit by GP98 of data by Ramaker & Peek (1976) is 
shown in Fig. [31 

2.4. Charge transfer, collisional excitation, and 
dissociation 

The H+-H2 and the H-H^ systems represent one of 
the most important collision complexes in the hydrogen 
plasma, called in short the H;^ collision system. In addi- 
tion to elastic collisions H2-H"'" and Hj-H, the relevant 
reactions are: (i) charge transfer. 



H+ 
-H 



H2(f) 

and (a) collisional dissociation. 



H+(i/)+H 



H^ 



f H- 
H+ 



(«), 
(6). 



The most complete theoretical work describing these 
chemical reaction pathways has been recently developed 
by Krstic and collaborators (Krstic 2002, 2003, 2005; 
Krstic et al. 1999, 2002) within the infinite order sudden 
approximation method (lOSA)^. As described by Krstic 
et al. (2002), charge transfer (a) is endoergic for v < 3, 
although with low threshold energies, and is strictly com- 
petitive with the vibrational excitation channel, espe- 
cially for those states than can overcome the kinetic bar- 
rier {v > 4). For higher vibrational levels, the process 

^ Data in electronic form are available at 
Ihttp: //www-cfadc.phy. ornl .gov/astro/ps/dat a/home .html I 
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Fig. 4. — H2-H"'" charge transfer rate coefficient. Triangles, total 
rate in LTE, summed over all possible initial and final states, from 
the cross sections of Krstic (2002) and Krstic et al. (2002, 2003); 
dashed curve, fit by Savin ct al. (2004) for }l2{v = 0); solid curve, 
fit listed in Table E] 
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Fig. 5.— H^-H (filled triangles) and H2-H+ ( open triangles) 
total dissociation rate coefficients in LTE, summed over all possible 
initial and final states, from the cross sections of Krstic (2002) and 
Krstic et al. (2002, 2003). The fits listed in Tableware also shown 
as solid and dashed lines, respectively. 

is exoergic, as the process of eharge transfer (6). for all 
V. The LTE rate coefhcient for the charge exchange re- 
action (a) is shown in Fig. |4l and compared with the fit 
by Savin et al. (2004) for H2(w = 0). For the reaction 
(6), cross-sections have been extrapolated at low ener- 
gies using a Langevin-type power law; the resulting rate 
coefficients have been scaled in order to obtain an LTE 
value corresponding to the experimental measurement by 
Karpas ct al. 1979 at Tgas = 300 K. The rate coefficients 

for the coUisional dissociation reactions (a) and (b) are 
shown in Fig. [5] The fits arc given in Table [3l 

2.5. H2/H vibrational-translational transfer 

Elastic collisions among hydrogen atoms and molecules 
deeply modify the initial vibrational distribution func- 
tion of molecular components in a hydrogen plasma. 
They are generally classified as V-T (vibrational- 
translational) and V-V (vibrational-vibrational) pro- 
cesses. In the former case, energy transfer between the 
translational degree of freedom of the colliding atom and 
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T[K] 

Fig. 6. — H2 dissociation by H impact: Triangles: rate coefficient 
computed from the cross sections of Esposito et al. (1999); solid 
line: fit reported in Table [3] 

the vibrational state of the target molecule occurs; the 
latter represents the energy exchange between the vibra- 
tional manifolds of the colliding molecules. In both cases, 
multi-quantum transitions can occur, coupling the entire 
set of vibrational levels. Recently, Krstic et al. (2002) 
have calculated V-T cross sections for collisions between 
H and H2, 

H + H2(i') ^H + H2(f'), 

using the quantum mechanical lOSA method. However, 
while data are available for the collision system 11+- 
Il2(v), no state-to-state information on collisions among 
hydrogen atoms and molecules is available for all {v,v') 
pairs. The V-T rate coefficients used in the present work 
have been obtained from cross sections computed by Es- 
posito et al. (1999, 2001) using the quasiclassical trajec- 
tory (QCT) method. 

Collisions of II2 with He atoms have been studied by 
e.g. Clark (1997), Flower et al. (1998), Balkrishnan et 
al. (1999), Lee et al. (2005), but no data for each {v,v') 
pair are available. In this work we have neglected V- 
T processes with He. This assumption is based on the 
observation that vibrating nuclei in H2 are lighter than 
the colliding He atom; moreover, the relative abundance 
of He is small: for these reasons, the V-T transfer with 
He can be considered less efficient than V-T transfer with 
H. Also molecule-molecule collisions have been neglected, 
owing to the small fractional abundance of H2 with re- 
spect to H. 

2.6. CoUisional dissociation of H2 

The cross sections for the dissociation of H2 induced 
by collisions with H, 

H-hH2(w) ^H + H + H, 

have been computed by Esposito et al. (1999, 2001) with 
the QCT method. The whole vibrational manifold of 
H2 has been explored. The total LTE rate coefficient is 
shown in Fig. [51 and the fit is given in Table [31 

2.7. Photoionization of H2 

Ionization of H2 can proceed through collisions with 
particles and photons. Data on ionization caused 
by electron impact are described by Liu & Sheman- 
sky (2004), and compared with photoionization data. 
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Fig. 7. — H2 photoionization rate coefficient. Triangles, total 
rate in LTE, summed over all possible initial and final states, from 
the cross sections by Flannery et al. (1977); dashed curve, fit by 
GP98 based on data by ONeil & Reinliardt (1978); solid curve, fit 
given in Table [3] 

Because threshold energies for the electron-impact ion- 
ization cross sections arc much higher than thermal en- 
ergies available at the redshift range considered in the 
present work (i? > 16 cV vs. E < 2 cV), we have con- 
sidered only photoionization of H2 , 

Several theoretical calculations of H2 photoionization 
have been proposed during the years (Lewis Ford et 
al. 1975, ONeil & Reinhardt 1978). This chemical re- 
action couples each vibrational level of H2 with each vi- 
brational level of . Available data on vibrationally re- 
solved {v, v') cross sections have been published by Flan- 
nery et al. (1977), following the theoretical and compu- 
tational procedures described by Tai & Flannery (1977). 
Fig. [7| shows the comparison between the total photoion- 
ization rate coefficient calculated in the hypothesis of 
LTE using the cross sections by Flannery et al. (1977), 
whose fit is given in Table [31 and the fit by GP98 based 
on data by ONeil & Reinhardt (1978). 

2.8. Photodissociation of H2 
Photodissociation of H2, 

can occur both directly (absorption from a lower level 
into the continuum) or indirectly (the so-called "Solomon 
process" ) , a process consisting in the absortpion into an 
individual level of a bound upper state, followed by dis- 
sociation. In the present work, we considered both de- 
struction pathways. For the direct process, vibrationally 
resolved rate coefficients are given by 

fcf (rrad)=4^ / ^\^MTr,d)d,^, (3) 
Jo '^^ 

where the cross sections a^^' are taken from the work 
by Allison & Dalgarno (1969). More recently, new 
calculations of photodissociation cross sections have 
been published by Glass-Maujean (1986) and Zucker & 
Eyler (1986), including excited electronic states; how- 
ever, no significant changes have been introduced for 
the Lyman and Werner systems. A thorough analysis 
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Fig. 8. — H2 indirect photodissociation rate coefficient. Tri- 
angles, total rate in LTE, both considering the Lyman and the 
Werner system, obtained using Eq. ^ and data by Dalgarno & 
Stephens (1970) and Allison & Dalgarno (1970); dashed curve, fit 
by S08 on Glover & Jappsen (2007) data; solid curve, fit given in 
Table [3] 

of the direct process with related fits and reaction rates 
for the kinetics of the early Universe is given in Coppola 
et al. (2011). 

For the calculation of the rate coefficient of the indirect 
process, we have used the discrete version of Eq. ([3]) 

fcJ?^'(T'rad) = ^A^„/„„e„ Ja„„ (Trad), (4) 

mcC^ 

where u is the vibrational level of both Lyman 
and Werner electronically excited intermediate systems, 
whose dissociation efficiencies have been taken from 
Dalgarno & Stephens (1970). The oscillator strengths 
fuv for the same molecular states have been calculated 
by Allison & Dalgarno (1970), and the vibrational en- 
ergies of both Lyman and Werner states by Fantz & 
Wiinderlich (2006). In Fig. |8] we show the comparison 
between the fit of the rate coefficient by S08 of the data 
by Glover & Jappsen (2007) and the total LTE rate co- 
efficient obtained using Eq. (|4]). The present calculation 
includes both Lyman and Werner systems, and the con- 
tribution of the total vibrational manifold has been con- 
sidered. 

2.9. Dissociative excitation and dissociative 
recombination of 

The fate of the dissociative scattering processes be- 
tween electrons and strictly depends on the energy 
of the collision. The effectiveness of these chemical re- 
actions is connected to the energy thresholds of the vi- 
brational states of the intermediate excited states: at 
high collision energies, the favored dissociative channel 
is dissociative excitation (DE): 

However, experimental evidence (Yousif & Mitchell 1995) 
has shown that large DE cross sections are operative at 
energies ^ 0.01 eV. For this reason, other reaction path- 
ways become important at these energies, involving elec- 
tron capture into excited electronic states of H2 . In par- 
ticular, H2(^S^) or auto-ionizing dissociative Rydberg 
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states lying below the Hj(2pcr„), both auto-ionizing in 
the continuum of H^, have been considered in the dy- 
namical mechanisms of dissociative excitation. 

At low energies, DE processes are competitive with 
dissociative recombination (DR), in which the incident 
electron is temporally captured into the Rydberg state 
by transferring its kinetic energy to vibrational motion. 
The unstable molecule can dissociate or go back to the 
ionizing state by an elastic or vibrational transition. This 
mechanism is sometimes called "indirect process" of DR. 
In the direct process, excited hydrogen molecules disso- 
ciate in the vibrational continuum of H2. In this case, 
two H atoms are produced, one in the ground electronic 
state, the other in states with principal quantum number 
n>2: 

Ilt{v)+e- ^ H + H(n > 2). 

Different studies have focused on DE and/or DR of 
Rl (e.g., Takagi 2002; Stroe & Fifirig 2009; Fifirig & 
Stroe 2008; Motapon et al. 2008; Takagi et al. 2009). We 
have used data from the work by Takagi (2002), where 
DR and DE cross sections have been calculated using 
the multi-quantum defect theory (MQDT). For each vi- 
brational level, Takagi (2002) provided both DR and DE 
cross sections. All excited dissociative states are included 
(all Rydberg states of 5 different symmetries) and for the 
lowest dissociative state, the electronic interaction is fully 
taken into account. Fig.[9]sliows the comparison between 
the total dissociative recombination rate coefficient cal- 
culated in the hypothesis of LTE using the cross sections 
by Takagi (2002), and the fit by GP98 based on data by 
Schneider et al. (1994). Our fit is given in Table El Al- 
though the cross sections of this process have a threshold 
component which corresponds to a repulsive channel, the 
threshold feature is located at too high energy to affect 
results in the present application. The main features rel- 
evant here are instead the presence of many resonances 
at low energy which convolute into a strongly decreasing 
trend for energies below a few eV, and a global increase 
of this low energy portion of the cross sections with v. 
As a result, the LTE rate coefficient decreases with Tgas 
until a few thousand K, where the raising population of 
vibrational levels reverses this trend. The fit by GP98, 
reported in Fig. [9] for reference, was based on data by 
Schneider et al. (1994), given up to Tgas = 4000 K and 
therefore showed a monotonic trend. 

2.10. Dissociative attachment 

Electron-impact inelastic processes of vibrationally ex- 
cited H2 molecules, 

ll2{v)+e- ^H + H-, 

play an important role in the kinetics of a low- 
temperature hydrogen plasma. Especially for inelastic 
collisions involving an electronic transition, two impor- 
tant features can be noted: if a vibrationally excited 
molecule is involved, the threshold decreases and the 
cross section of the process increases (for some processes 
dramatically) with increasing vibrational excitation of 
the molecule. In this work we have used data calculated 
by Celiberto et al. (2001), based on an improved version 
of the resonant scattering model originally developed by 
Fano (1961). The dynamics of this route provides for 
the formation of a temporary negative molecular ion. 
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Fig. 



9. — dissociative recombination rate coefficient. Tri- 



angles, total rate in LTE, summed over all initial states, from the 
cross sections by Takagi (2002); dashed curve, fit by GP98 based 
on data by Schneider et al. (1994), given in the temperature range 
20 K< Tgas < 4000 K; solid curve, fit reported in Tablc[3] 



whose fate is either ionization or dissociation. The fit 
of the LTE rate coefficient can be found in Capitelli ct 
al. (2007). 

2.11. Electron collisional excitation 0/H2 

The process of vibrational excitation of molecules by 
electron collisions, 

H2(w) + e" ^ H2(w') + + hv, 

also known as the E-V process, selectively populates high 
vibrational levels and represents one of the main non- 
equilibrium pathways in the vibrational kinetics of hy- 
drogen plasma. It is a two-step process that links a vi- 
brational state of the electronic ground state to another 
vibrational state of the same manifold, via an interme- 
diate singlet state that can radiatively decay. As in the 
case of dissociative attachment, we have used data from 
Celiberto et al. (2001). 

2.12. Formation o/HJ via HeH^ 
The HeII+ channel for the formation of H J [v) , 

HeH+ + H ^ H+(i;) +He, 

has been considered in several studies (e.g. Dal- 
garno 2005, HP06, S08). The rate coefficient usually 
adopted in these kinetics models is estimated from the 
experimental work by Karpas et al. (1979) or data by 
Linder et al. (1995). However, no vibrationally resolved 
information is available. In order to overcome this limita- 
tion, we follow here an approach similar to that of Hirata 
& Padmanabhan (2006), where this reaction is assumed 
to equally populate all levels compatible with energy bal- 
ance, i.e. all levels with energy E < —1.844 eV with re- 
spect to the dissociation limits. Since in the present 
work individual rovibrational levels are not considered, 
wc follow a simplified approach by assuming that this re- 
action channel produces molecules in the first three 
vibrational levels with equal rates. 

2.13. Radiative transitions 0/ II2 



H2 and in the primordial gas 



7 



Quadrupolar transition probabilities for H2(v) have 
been calculated by Wolniewicz et al. (1998)''. The 
authors improved earlier calculations by Turner et 
al. (1977) by adopting a more accurate potential energy 
curve, providing tables for the spontaneous decay from 
an initial to a final vibrational level for A J = 0, ±2. The 
rovibrational energy levels used to perform the sum of 
available transition rates have been calculated using the 
WKB method (Esposito 2010, priv. comm.). They are in 
very good agreement with classical calculations by Kolos 
& Wolniewicz (1964). 

2.14. Radiative transitions of 

Calculations of quadrupole transition probabilities for 
HJ(?;) follow the method used by Wolniewicz et al. (1998) 
for H2. Radiative transitions for have been calcu- 
lated by Posen et al. (1983) for vibrational-vibrational 
transitions at given J and A J. The energies of the rovi- 
brational levels used to perform the sum have been cal- 
culated by Hunter et al. (1974). 

3. CHEMICAL NETWORK 

To describe the time evolution of the chemical species 
during the expansion of the Universe, we solve numeri- 
cally the system of ODEs 



driy 



(5) 



where Ry^i are the spontaneous and stimulated radiative 
rate coefficients; Pyyi are the destruction rate coefficients 
of the w-th species by photons; C^^i are the destruction 
rate coefficients for the u-th species for collisions with the 
w'-th chemical partner; and are the formation rate 

coefficients of the u-th species due to collisions between 
the w'-th and v"-Xh species, photons included. 

The radiative rate coefficients are related to the Ein- 
stein coefficients by 



Rvv' 



Trad) iiv'<V 
Tiad) if v' > V, 

^ ( Ayy' [1 + r]{u„y> , T^ad)] if v' < V 

\ 5^'A„/„?7(M„„/,T,.ad)/.9t, a v' > V 

where Vw' is the frequency of the transition u — v', 
^t, = 1 is statistical weight of the w-th level, u{i/y 

v' 1 Tj-ad) 

is the Planck photon distribution, and ?y(^'i,t)' , Trad) — 
[exp(/i^'^„//fcTrad) - 1]"^ 

The chemical network is completed by the equations 
for the gas and radiation temperature evolutions and the 
equation for the redshift. 



dt 
dz 



1 



{i + z)H{zy 



(6) 



where 



H{z) = Hoy^Q,{l + z)4 + a,i(l + z)^ + ^k{1 + 2)2 + f^A. 

(7) 



^ Data in electronic format are ava laible 
|http: //cfa-MW«. harvard ■ 6du/~ simbotin/4pole .html[ 



The number density of hydrogen atoms is 
3Hi 



-(l-rp)(l + z)^ 



(8) 



SttCtoh 

where f2(, is the baryon fraction, Yp is the helium mass 
fraction, G is the constant of gravitation and mn the 
atomic hydrogen mass. 

Numerical values of the cosmological parameters Hq, 
Tq, fir, ^m, ^b, and Yp have been obtained from 
WMAP5 data (Komatsu et al. 2009, see Table [2]). The 
initial fractional abundances for H, He and D are also 
listed in the same table. The hydrogen ionization frac- 
tion was computed using the routine RECFAST (Seager 
et al. 1999, Seager et al. 2000, Wong et al. 2008). Hehum 
was assumed to be fully neutral in the redshift range of 
interest (2 < 2000). The electron density was determined 
by imposing the condition of charge neutrality. Atomic 
and molecular weights were taken from the NIST web- 
page ®. 



TABLE 2 
Cosmological parameters 



Cosmological parameter 


Numerical value 




100 ft km s-i Mpc-i 


h 


0.705 




3141 


To 


2.725 K 




0.228 


Hb 


0.0456 








f^m/(l + Zcq) 


Ha 


0.726 




1 — — ^^m ~ ^^A 


yp 


0.24 


/h 


0.924 


/hc 


0.076 


/d 


2.38 X 10"^ 



at 



Data tables of vibrationally resolved rate coefficients 
have been inserted as input of the routine LSODE ^ 
used to integrate the chemical system, which consists of 
47 differential equations, one for each chemical species 
introduced in the model and one for the gas tempera- 
ture. An implicit method has been used to perform the 
integration, being the kinetic problem stiff. Linear in- 
terpolation of the rate coefficients and of the ionization 
fraction is performed in logarithmic scale at each step of 
integration. 

4. RESULTS AND DISCUSSION 

The resulting fractional abundances of H2 and for 
each vibrational level are shown as a function of redshift 
in Fig. [TUl The total fractional abundance of H2 shows 
a rapid increase at three epochs: at z « 1500 by the 
charge transfer channel, dominant at high temperatures, 
at z ?a 300 by radiative association formation, that 
modulates the charge transfer channel, and at z w 100 
by the associative detachment process. At z = 10 the 
fractional abundances of H2 and are 5.76 x 10~^ and 
6.56 X 10~^^, respectively. Our value of the H2 abun- 
dance is in good agreement with that obtained by HP06 

* |http : //webbook. nist .gov/ chemistry /name-ser. html I 

® https : / /computation. llnl .gov/casc/odepack/download/lsode.agree .htm] 
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Fig. 10. — Vibrational level populations of H2 {top panel, from 
V = to V = 14) and {bottom panel, from v = to v = 18) as 
function of redshift. In Both cases, the solid curve shows the total 
fractional abundance. 



at z 



20) but about twice the 
10). As for H+, our 



(/(H2) 6.0 X lO-'^ at z 
value of V09 (2.7 x lO"'^ 
abundance is in good agreement with the result of V09 
(/(H+) = 6.7 X 10"" at z= 10), whereas HP06 obtain 
three different abundances depending on the value of the 
H^-H charge exchange reaction. 

For H2, a marked non-equilibrium distribution of pop- 
ulations is established at z « 300, followed by a plateau 
involving levels from v = 1 to 9. Levels above v = 9 are 
excluded from this last plateau as their populations drop 
dramatically due to the endoergic character of the pro- 
cess of associative detachment (see Section [^?T|) . There- 
fore, when the thermal kinetic energy is much lower than 
the energy thresholds, the formation of highly excited 
{v > 10) H2 molecules is strongly suppressed. Such a 
strong splitting of level histories is not observed at higher 
z, since the main H2 formation mechanism (H + 
charge transfer) is exoergic for all v (in constrast with the 
charge transfer H+ + H2 , that is a threshold process) . 

The H2 vibrational distribution functions obtained at 
different values of the redshift are shown in Fig. 111! com- 
pared with the equilibrium curves corresponding to the 
Boltzmann vibrational distributions at the correspond- 
ing value of the gas temperatures. The shape of the level 
population distribution of H2 at low z can be understood 
and even roughly evaluated on the basis of a balance be- 
tween radiative and formation rates. Indeed, the bunch- 
ing of population of excited levels for low z is due to the 
fact that each population is determined essentially from 
the ratio of associative detachment and radiative rates. 
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Fig. 11. — Vibrational level populations of H2, normalized to the 
total H2 fractional abundance, at various redshifts. The Boltzmann 
distributions at each redshift are indicated with the label "B" . 

The ratios for all exothermic channels {v <10) have com- 
parable values on a enlarged logarithmic scale like that 
used here. 

For HJ, the most relevant formation process is ra- 
diative association. The vibrational resolution of this 
species shows that the formation process, being faster for 
highly excited vibrational levels of the products, leads to 
a pronounced tail in the vibrational distribution (Fig. 
[T2|) . The latter is characterized by a high population of 
the first few vibrational levels followed by a suprathermal 
tail, in qualitative agreement with the results of HP06 
(see their Fig. 4). However, for HP06 the dominant for- 
mation process of H J is via HeH+ . This result supports 
the claim of HP06 that the radiative association channel 
plays a minor role among formation pathways in their 
model. Also in this case, we observe a strong group- 
ing of the fractional abundance of excited levels which 
is explained by the balance between state specific level 
production rates and radiative transitions to w = 0. 

These findings suggest that the processes leading to 
redistribution of the vibrational quanta need further 
attention. This study is also appropriate here be- 
cause such processes cannot be modeled in the frame- 
work of the usual, not state-resolved approach. In our 
case, such redistributing channels are V-T and sponta- 
neous/stimulated radiative processes. In order to un- 
derstand how these pathways effect the vibrational dis- 
tribution of H2 and HJ, further numerical experiments 
have been performed in various regimes, corresponding to 
particular physical conditions. Fig. [T^ shows the results 
obtained for the case in which radiative processes have 
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Fig. 12.— Same as Fig. [ll] for H^. 

been omitted. The figure shows that, aUhough radiative 
processes are not fast enough to thermahze level popula- 
tions, they are essential to produce the actual vibrational 
distribution, and their neglect leads to complete different 
results especially at low z. 

On the other hand, removal of V-T processes does not 
produce any appreciable variation with respect to the 
results of Fig. [TOl since the rates of V-T processes are 
much smaller than the radiative ones. The role of V- 
T processes is better appreciated by looking at Fig. UM 
where V-T processes and radiative decay are both ig- 
nored. Some variations with respect to Fig. [13] are ob- 
served: a higher vibrational temperature, a different high 
V plateau for H2, especially at low z, and some differences 
of the total fraction. 

The main destruction channel (dissociative recombi- 
nation) active at lower z is deeply affected by radiative 
processes: when these are ignored, the destruction chan- 
nel becomes more efficient and the fractional abundance 
at the second peak (at z ^ 300) is reduced by about a 
factor of 2. As a further check, removing radiative pro- 
cesses together with the loss channel leads to small vari- 
ations of the fractional abundance. Thus, redistribution 
among levels is able to affect even the results for the to- 
tal H^, confirming the relevance of the non-equilibrium 
vibrational distributions. 

No rotationally resolved data exists for most of the for- 
mation and destruction processes included in the present 
work, so a state-to-statc kinetics cannot at the present 
stage go beyond the resolution of molecular vibrations. 
For H2 and , the typical rotational energies are of the 
order of 0.01 eV: this suggests that the way in which ro- 




FlG. 13. — Same as Fig. 1101 including formation and destruction 
channels, V-T processes but no radiative transition. 

tation is included in the kinetic model can affect the low 
z trend of molecular abundance. To test the sensitiv- 
ity of the chemical abundances to the rotational struc- 
ture of H2 and H^, we have computed the radia- 
tive association rate coefficient (Subsection 12. 3p assum- 
ing that the rotational manifold for each vibrational level 
reduces to only one J, cither the lowest or the highest. 
Even if this hypothesis implies a very strong rotational 
non-equilibrium, it has negligible effect on the freeze-out 
value of H^, although the abundance of this species at 
z ~ 1000 is changed. Furthermore, since the photodisso- 
ciation cross section is not rotationally resolved, the ro- 
tational degree of freedom only enters in Eq. ([2]) through 
Zrot- As a consequence, the nascent vibrational distri- 
bution function, which is the most important quantity 
here, is only marginally modified. 

5. CONCLUSIONS 

In this work we have performed calculations of the vi- 
brational distribution of both H2 and 112" for the con- 
ditions expected in the early Universe and based on a 
comprehensive state-to-state chemical kinetics. The vi- 
brational level distribution for these two species arc re- 
ported for a wide and continuous range of the redshift 
parameter. The results can be summarized as follows: 

1. The vibrational distribution function of H2 and 
assumes a quasi-equilibrium shape at redshift 
z ^ 1500; after that, extended plateau in the vi- 
brational level distributions form, underlying the 
presence of pumping phenomena for the interme- 
diate vibrational levels; full thcrmalization is not 
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Fig. 14. — Same as Fig. 1101 including formation and destruction 
channels, but no V-T process and radiative transition. 



observed, because vibrational relaxation processes 
are not fast enough to balance the strong vibra- 
tional selectivity of formation rates. 



2. Radiative processes play a fundamental role in the 
redistribution of vibrational quanta, affecting the 
vibrational distribution function and the overall 
fractional abundance (as in the case of H^); 



3. All these features can not be described in terms of 
the LTE distributions which arc usually assumed 
in chemical networks for the primordial Universe. 
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knowledge Jonathan Tennyson and loan F. Schneider for 
useful discussions on electron— molecule collisions, and 
MIUR and Univcrsita dcgli Studi di Bari, that partially 
supported this work ("fondi di Ateneo 2010"). 



APPENDIX 

ANALYTICAL FITS OF LTE RATE COEFFICIENTS 

In the present section, the analytic expressions of the vibrational LTE rate coefficients of the chemical processes 
introduced in the present model are reported. For processes that are vibrationally resolved in the final states, the sum 
on the entire final manifold has also been carried out. 
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TABLE 3 

Chemical reactions and their LTE rate coefficients 



Chemical process Rate eoefficient (m"* s ^) or (s 



1] H + n- -^H2(v)+e- logio fc = En=0 «2„ [logio (Tga. / 10^ )] 2" 

ao = —14.4 
a2 = -0.15 
ai = -7.9 X 10-3 

2] H2(t;)+H+->H+(t)') + H Infc = ao +aiTgas +a2/Tgas +a3Tg\,, 

ao = -33.081 
ai = 6.3173 X 10"^ 
a2 = -2.3478 x 10* 
as = -1.8691 X 10^3 

3] H+(j;) + H^H++H + H Infc = ao + aiTgas + a2/T'gas + asTg^^, 

ao = -32.912 
ai = 6.9498 x IQ-^ 
02 = -3.3248 X 10* 
as = -4.08 X 10-9 

4] H2(i;) + H+ ^ H + H + H+ Infc = ao + aiTgas + a2/Tgaa + a3T|^,, 

ao = -33.404 
ai = 2.0148 X 10"* 
a2 = -5.2674 x 10* 
as = -1.0196 X 10"* 

5] H2(«) + H->H + H + H fc = 1.9535 X 10-"'Tg-°-332®''exp{-4.9743 X 10* /Tgas) 

6] ll2iv) + hiy ^U+{v') + e- fe = 3.06587 X 10^ exp(-l. 89481 X lO^/Tgas) 



7] H2(i') + /i!^ -> H + H (indirect) Infc = 17.555 + 7.2643 x lO"'' Tgaa - 1.4194 x lO^/Tgas 

8] H+(,;) + e-^2H ^ = E'=o «-7^g"as, 

ao = 4.2278 X 10"" 
ai = -2.3088 x 10"^^ 
a2 = 7.3428 x IQ-^i 
as = -7.5474 x lO'^^ 
ai = 3.3468 x IQ-^s 
as = -5.528 X 10'=** 

9] H- + -> H + e- k = O.llT^J^ exp(-8.823 x 10^ /T,^d) 

fc = 8.0 X 10-*Tl |exp{-2.3 X lO^/T^ad) (non-thermal) 
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